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(~ | ■ We study the performance of Bacon-Shor codes, quantum subsystem codes which are well 

Q ' suited for applications to fault-tolerant quantum memory because the error syndrome 

can be extracted by performing two-qubit measurements. Assuming independent noise, 
we find the optimal block size in terms of the bit-flip error probability px and the 
phase error probability pz, and determine how the probability of a logical error depends 
on px and pz- We show that a single Bacon-Shor code block, used by itself without 
concatenation, can provide very effective protection against logical errors if the noise is 
highly biased (pz/px 2> 1) and the physical error rate pz is a few percent or below. We 
also derive an upper bound on the logical error rate for the case where the syndrome 
data is noisy. 

Keywords: Quantum error correction 

1 Introduction 

Bacon-Shor codes [TJ [5] are quantum subsystem codes [3J H] which are well suited for ap- 
plications to fault-tolerant quantum memory [SJ [3] , because error syndrome information can 
be extracted by measuring only two-qubit operators that are spatially local if the qubits are 
arranged in a two-dimensional lattice. In this paper we assess the performance of these codes. 

We consider noise models such that qubits in the code block are subject to both bit flip 
(X) errors and dephasing (Z) errors, where the bit flips occur with probability px and phase 
errors occur with probability pz- We assume that the noise acts independently on each qubit, 
and that the X and Z errors are uncorrelated. Under these assumptions we find the optimal 
block size of the code, and the failure probability achieved by this optimal code. We obtain 
analytic formulas for the optimal failure probability for the case of unbiased noise (px = Pz) 
and the case of highly biased noise (b = pz/px 3> 1). In both cases the formula applies in 
the asymptotic limit of small pz- Our results show that the failure probability of the optimal 
code falls exponentially in l/pz- We also consider the case where the code's syndrome bits 
are prone to error, and derive upper bounds on the failure probability in that case. 

We find that for the case of unbiased noise (pz = px = p), the optimal Bacon-Shor code 
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achieves the failure probability 



2 \ 1/2 /ln 2 2\ 1/0 / In 2 2 



BSOptimalFailProb(p) = ( " ) exp ( ) p 1/2 exp ( + 0(p) 



(1) 



or 



\a[BSOptimalFailProb{p)] = -A/p - (0.5) ln(l/p) + C* + 0(p), (2) 

where 

A = .0600566, C = .0175217. (3) 
For the case of highly biased noise, we find 

\n[BSOptimalFailProb(p z , b)] = ~A{b)/p z - (0.5) ln(l/p z ) + C(b) + 0{p z polylog(fc)), (4) 
where 

Mb) = l(w(Vb)) 2 + 0{b- 1 ' 2 \nb); (5) 
here W denotes the Lambert W function, with asymptotic expansion 

W(Vb) =lnVb-lnlnVb+ ^5 + O ( . (6) 

\nVb \ In 2 Vb J 

We also compute the asymptotic form of C(b) for b S> 1. 

In the case where the error syndrome is noisy, the reliability of the syndrome can be 
improved by measuring it repeatedly. We consider an idealized noise model such that qubit 
errors and syndrome measurement errors are equally likely and independent. For that model 
we derive a lower bound on the coefficient A{b) of 1/pz in the natural logarithm of the logical 
failure rate, finding 



where fi = 2.63816 is the connective constant of a self-avoiding walk on a two-dimensional 
square lattice, and Z(b) is a slowly varying monotonic function which ranges from 1/e to 1 
as b increases from 1 to infinity. This bound applies for any value of b > 1. 

Previous work p] |S] on the applications of Bacon-Shor codes to fault-tolerant quantum 
computing has focused on relatively small codes used at the bottom layer of a concatenated 
coding scheme. Here we emphasize that a sufficiently large Bacon-Shor code, used by itself 
without concatenation, can also be quite effective if the physical error rate is low enough. For 
example, if the syndrome is perfect, then the probability of a logical failure is below 2 x 10 -19 
for pz = .01, b = 100, and below 10~ 12 for pz = .03, b — 1000. Fault-tolerant circuits based 
on Bacon-Shor codes will be more fully discussed and analyzed in a separate paper [5], where 
we consider in particular the case of highly biased noise. 
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2 2D Bacon-Shor code 

The two-dimensional Bacon-Shor code [TJ [2] is a quantum subsystem code (really a family 
of codes), which pieces together two dual quantum repetition codes, one protecting against 
bit flip errors and one protecting against phase errors. Like any subsystem code [31 0], it can 
be defined by its "gauge algebra" — a set of Pauli operators that commute with the logical 
operators. The center of the gauge algebra is the code's stabilizer group, and the code space 
is determined by fixing the eigenvalues of all the stabilizer operators to be (say) +1. Other 
operators in the gauge algebra act nontrivially on "gauge qubits" but trivially on the code's 
protected qubits. 

Specifically, we consider combining a length-m repetition code to protect against Z errors 
with a length-n repetition code to protect against X errors, where Z and X denote the single- 
qubit Pauli operators az and ax- The resulting Bacon-Shor code block contains mn qubits 
arranged at the sites of an m x n square lattice. The gauge algebra is generated by two-qubit 
operators XX acting on pairs of neighboring qubits in the same horizontal row of the lattice 
and two-qubit operators ZZ acting on neighboring qubits in the same vertical column. There 
is one protected qubit, and we may choose the logical Pauli operators to be X = X® 11 acting 
on all the qubits in a column and Z = Z® m acting on all the qubits in a row. The stabilizer 
group has m + n — 2 generators — there are m — 1 X-type generators, each a product of 2n 
Xs acting on all the qubits in a pair of neighboring columns, and n — 1 Z-type generators, 
each a product of 2m Zs acting on all the qubits in a pair of neighboring rows. (We say that 
a Pauli operator is X-type if it is a tensor product of Xs and identity operators, and we say 
it is Z-type if it is a tensor product of Zs and identity operators.) Aside from the protected 
logical system, there are (m — l)(n — 1) gauge qubits. 



x x 

Fig. 1. Bacon-Shor code block for m = 5 and n = 3. A weight-(2m) Z-type stabilizer generator is 
the product of m weight-two Z-type gauge operators (solid red); a weight-(2n) X-type stabilizer 
generator is the product of n weight-two X-type gauge operators (hashed blue). 



Though each Z-type stabilizer generator (or "check operators" ) has weight 2m, its value 
can be ascertained by measuring m weight-two gauge operators. For example, the product 
of 2m Zs in the first two rows is equivalent to the product of m Z-type gauge operators, 
each acting on a pair of qubits in the first two positions of the same column. To determine 
the check operator, we may measure these m gauge operators (obtaining either +1 or —1 for 
each outcome) and multiply the results. Similarly, the value of a weight-2n X-type check 
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operator can be found by measuring n weight-two X-type gauge operators and multiplying 
the results. Furthermore, since all check operators commute with all gauge-qubit operators, 
measuring X-type gauge qubit operators to determine the value of an X-type check operator 
does not disturb the values of any Z-type check operators (though it may flip the value of 
Z-type gauge operators), and vice- versa. 

If m and n are both odd, then the code can correct (m — l)/2 Z errors and (n — l)/2 X 
errors. To describe the error recovery procedure, it is useful to exploit the freedom to "fix the 
gauge" , i.e. apply gauge operators to the block that commute with the logical qubit operators. 
An even number of X errors in a row can be "gauged away" by applying gauge operators. If 
there are an odd number of X errors in a row, we may "gauge away" all errors except one, 
which can be moved into the first position in the row. Thus, in an appropriate gauge, X 
errors occur only in the first column. Likewise, in an appropriate gauge, Z errors occur only 
in the first row. In this gauge, then, the check operators test whether two neighboring qubits 
in the first column agree in the Z basis, and whether two neighboring qubits in the first row 
agree in the X basis. Error recovery proceeds by applying X to at most (n — l)/2 qubits in 
the first column, to restore all Z-type check operators to the value +1, and by applying Z to 
at most (m — l)/2 qubits in the first row, to restore all X-type check operators to the value 
+1. Thus, error recovery is successful if the number of rows with an odd number of X errors 
is at most (n — l)/2 and if the number of columns with an odd number of Z errors is at most 
(m- l)/2. 

If we choose the gauge so that ZZ = +1 for any pair of qubits in the same column, then 
the eigenstates of X with eigenvalues ±1 become 

|±> c cx(|00---0)±|ll...l}f m ; (8) 

these are tensor products of m lcngth-n "cat states" in the standard basis, one for each 
column. If we choose the gauge so that XX = +1 for any pair of qubits in the same row, 
then the eigenstates of Z with eigenvalues ±1 becomes 

|0)c cx (|+ + ...+) + | )f\ 

|l)c cx (|+ + . .. + )-| )f n - (9) 

these are tensor products of n length-m "cat states" in the dual basis, one for each row. 

To perform a destructive measurement of the logical operator Z , we can measure all mn 
qubits in the Z basis, compute the parity of the outcomes for each row, and then decode 
the measurement by performing a majority vote on the row parities. Likewise, to perform a 
destructive measurement of the logical operator X, we can measure all mn qubits in the X 
basis, compute the parity of the outcomes for each column, and then decode the measurement 
by performing a majority vote on the column parities. 

3 Failure probability for the 2D Bacon-Shor code 

Consider an independent noise model in which the probability of an X error for each qubit 
is px, and the probability of a Z error is pz (that is, for each qubit, there is no error with 
probability (1— px)(l — Pz), an X error with probability px{^—pz), a Z error with probability 
(1 —px)pz and an iY = ZX error with probability pxPz)- For any fixed values of px and pz, 
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there is an optimal choice for the code dimensions to and n which minimizes the probability of 
a logical error. We will estimate the optimal failure probability in the limit of small p. Since 
X and Z error recovery are identical (except for the interchange of the rows and columns of 
the lattice), we will estimate the probability of a logical Z error; the same derivation also 
applies to the probability of a logical X error. 

First we note that for a column of length n, if Z errors occur independently with probability 
p at each position, then the probability of an odd number of errors in the column is 

OddProb{ Pl n) = ^ (1 - (1 - 2p) n ) . (10) 

We derive this formula by observing that the terms of even order in p in the binomial expan- 
sions of ((1 — p) + p) and ((1 — p) — p) are identical, while the terms of odd order in p are the 
same except for a sign flip; hence the difference | (1 — (1 — 2p) n ) sums up all the odd-order 
terms. To minimize the probability of a logical error for given p, we will choose n and to so 
that p 2 mn is 0(1). Therefore in the limit of large n we may use the approximation 

OddProb{p, n) = ^ (l - exp [-2p(l + p)n + 0(p 3 n)] ) . (11) 

We keep track of the 0(p 2 n) in the exponential because OddProb(p,n) will be raised to a 
power 0(m) in the expression for the Bacon-Shor failure probability. 

Now consider a length-m repetition code, where bit errors occur independently with prob- 
ability x. If m is odd, an encoded error occurs if the number of bit errors is (to + l)/2 or 
more. Therefore the failure probability of the repetition code is 

RepFailProb(x,m) = ^ ( ^ j x k (l - x) 

fc=(m+l)/2 



1-x 

From the Stirling approximation, 



k 

1/2 (m-l)/2 



m \ _ / 2 \ 1/2 / 2 



m+l 



neglecting a multiplicative 0(l/m) correction. Making another 0(1/to) multiplicative error, 
we may replace the exponential inside the sum over r by 1, obtaining 

/ o \ 1/2 / \ 1/2 (m-l)/2 , . r 

RepFailProb(x,rn)«2-^—j [x(l x)]^ ^ \JL-} , (14) 



and we also make a negligible error by extending the upper limit on the sum to infinity, finding 
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and thus 



RepFailProb(x 



To compute the probability of a Z-type logical error for the Bacon-Shor code, we substitute 
OddProb(p,n) for x in the expression for RepFailProb(x , to) , finding 



[4x(l - x)] m/2 f 

and 
hence 

BSZFailProb(pz, m, n) 



(l - cxp [-4p(l +p)n + 0(p 3 n)])' 



exp 



2p 2 mn[l + O(p) + 0(p 2 
e ipn - 1 



z/2 



(1 



4pny»/2 



(1 - 2z) 2 ~ 4 1 j ' 



RepFailProb(OddProb(pz , n), to) 



(17) 



(18) 



1 



V2 



7TTO 



: exp 



2p 2 z mn 



(e 4pzn - 1) V2 (1 - e - ipzn y 



i/2 



(19) 



up to multiplicative corrections higher order in pz (assuming pzn and pzm are 0(1)). By 
the same argument, the probability of an JT-type logical error is given by a similar expression, 
but with m and n interchanged and pz replaced by px- 



BSXFailProb(p x ,m,n) = RepFailProb(OddProb(p x ,m),n) 



'2im 



cxp 



2p x mn 



n/2 



(20) 



up to multiplicative corrections higher order in px (assuming pxni and pxn are 0(1)). 
3.1 Unbiased noise 

If the noise is unbiased (px = pz = p), then the optimal Bacon-Shor code is symmetric 
(m = n), and the X-type and Z-type logical errors occur with equal probability. To find the 
optimal value of n for small p, we note that 



In [BSFailProb{p, n)] « | In (l - e - Apn ) + 



(21) 



where the ellipsis indicates corrections suppressed by powers of p or 1/n. This expression, 
regarded as a function of n, attains its minimum when y — pn satisfies 



In (1 - e- Ay ) + 



4y 



e^v - 1 



= 0, 



(22) 



or Ay = In 2. (Though n is actually required to be an odd integer, we may ignore this 
requirement if p is small and n is correspondingly large.) Substituting Apn w In 2 into the 
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expression for the Bacon-Shor failure probability Eq.(|19j) we find 

/ 2 \ 1//2 /ln 2 2\ / In 2 2 

BSOptimalFailProbOp) ps I — — ) exp ) p 1/2 exp 

\7rln2/ \ 8 / 

fi mwart r- f -0600566 \ 
= (1.01768) exp , (23) 



ln[BSOptimalFailProb(p)} w -A/p + B Lap + C + 0(p) (24) 

where 

A = .0600566, B = .5, C = .0175217. (25) 
This optimal failure probability is achieved by choosing the linear size n of the code such that 

pn = i In 2 = .173287. (26) 

In terms of the asymptotically optimal linear size n = Jj- 2 - , the optimal failure probability can 
be expressed as 

, , , / 2 \ 1/2 /ln 2 2\ /ln2\ 1/2 ( nln2 

BbOptimalFailProbyn) ~ I — j — - I exp I — — J I — — I exp 

oln2/8 

= ?- n/2 . (27) 



V27rn 

Although there is strictly speaking no accuracy threshold for this family of codes, the 
performance of the optimal code is very good when the error rate is sufficiently small. For 
example, for p — .001, the optimal value of n is 173, and the corresponding failure probability 
is 2.638 x 10" 28 , while the asymptotic formula Eq.flU]) predicts 2.663 x 10~ 28 . 

3.2 Highly biased noise 

If the noise is highly biased (b = pz/px ^ 1), then the optimal Bacon-Shor code is highly 
asymmetric — it combines a length-m repetition code to protect against Z errors with a 
length-n repetition code to protect against X errors, where m 3> n. In the limit of large noise 
bias, we can find analytic formulas for the optimal values of m and n and for the corresponding 
optimal probability of a logical error. 
From Eq.finj) and Ea.([2"0"]). we have 

ln[BSZFailProb(pz,m,n)] y In (l - e~ 4pz ") + • • • , 
In [BSXFailProb(p x ,m,n)] w - In (l - e - 4pxm ) + ■ ■ ■ . (28) 

where the ellipsis indicates corrections suppressed by powers of pz, px, m -1 , For 
b = pz/px 3> 1) the optimal values of m and n are such that Apxm and e _4pz ™ are both 
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small; hence we can justify keeping the leading terms in a power series expansion in these 
small quantities, obtaining 

hilBSZFailProb^pzTm^n)] m — — cxp (— Apzn) — — exp (— 8pzn) + ■ ■ ■ , 
ln[BSXFailProb(px,m,n)] sa — hi (Apxrn) — pxTnn +■■■ . (29) 

For now we neglect the non-leading terms in the expansion of both logarithms, which we will 
verify a posteriori are additive corrections of order 6 -1 / 2 ln&. 

We will minimize the failure probability subject to the constraint that (the leading contri- 
butions to) the logs of the Z and X failure probabilities are equal. Because we are neglecting 
the higher order terms in both Eq. (|2"5| and Eq. (|2"9"l) . the solution we find many not be the true 
optimum, but we will see that it provides an upper bound on the optimal failure probability 
which is reasonably tight provided pzln b and 6 _1 / 2 ln6 are small. Defining the variables 

X = 4p x m, Z = 4p z n, (30) 

we want to find the values of X and Z that minimize the function 

F{X,Z) = Z\nX (31) 

(the log of the failure probability multiplied by 8pz — Sbpx) subject to the constraint 

Z\nX = -bXe~ z . (32) 

Introducing a Lagrange multiplier A to do the constrained minimization, we obtain the equa- 
tions 

-(\-l)Z = XbXe- z = (A- l)1nX (33) 
together with the constraint equation Eq. (132|) : these equations imply 

X = e~ z and Ze z = Vb => X = Z/Vb. (34) 
Evaluating F at its minimum yields 

i^min - -Z 2 . (35) 

The solution to the equation Ze z = Vb is the Lambert W function Z = W(y/b), which 
has the asymptotic expansion 

Z = W(Vb) = In Vb In In Vb + + O ( (36) 

In v b \ In v b J 

for b ^> 1 , and the log of the optimal failure probability is 

W 2 (Vb) 

In [BSZFailProb(p z , m, n)] w In [BSXFailProb{p z , m, n)] « - + ■■■. (37) 
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The optimal code has dimensions 

n = « -^—W{Vb) a -L hi A, 

m = - — « — VbW(Vb) sa — A In A (38) 
4px 4pz 4^2 

with aspect ratio m/n = A- As the bias increases with pz fixed, the code size creeps up 
slowly in the X-protection direction, and more rapidly in the Z-protection direction. 

For the optimal code the probability of an odd number of X errors in a row decreases as 
b increases according to 

OMProb^m) « \ (1 - ) « Px m « « 1^ « (39) 

while the probability of an odd number of Z errors in a column asymptotically approaches 
1/2 according to 

OddProb(p z , n) a ~ (l - e" 2 ^") = i (l - e^ 2 ) a i (l - y/x) , (40) 

or 

1 - OddProb(p z ,n) « 1 f . (41) 



2 ; 2 I A 



The optimal failure probability is much smaller for b 3> 1 than in the case b = 1 for the same 
value of pz, but the price paid is that the block size is correspondingly significantly larger: 



VbW 2 {Vb) 

mn = . 

16p| 

as compared to 



(42) 



" 2 = S <«> 

in the unbiased case. 

Using Eq. (p8|) . we can now evaluate the subleading terms in Eq. ([29)) . finding 



— exp {-8p z n) ps — -A 



2 



i /A w(A) \ /w(v^)\ _ vk 3 (A) 



4 4 4\ v 4p Zy /\ y A/ 16pzA 

-pxmn « — i — - — p — - = y —=. (44) 

b \ Ap z ) \ ±pz J IQpz^/b 

Aside from logarithmic factors, these terms are suppressed by 0(6 _1 / 2 ln&) compared to the 
leading terms in Eq. (|2"9")l . We therefore expect our estimate for the optimal code dimensions to 
be accurate for b 3> 1. However, although these corrections are small compared to the leading 
terms, they are not necessarily negligible, and in particular they can become important as 
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pz — > with b fixed. To obtain accurate results for small pz, we should resum the power 
series expansion in Eq.(|29[). or in other words use the full expressions in Eq. (|28|) . with the 
appropriate values of m, n, Z, and X plugged in. Thus, 



BSZOptimalFailProb(p z , b) 



(l _eY«(l-M) 



VbW(Vb) 
8PZ 



Vb 



BSXOptimalFailProb(p z , b) (l - e" x )™ /2 « (l - cxp ( 



_ p -w(Vb) 



W(Vb) 



, (45) 



provide more accurate approximations than exp ^— W 2 (Vb)/8pz 

Plugging our solutions for Z and X back into the expressions for the prefactors in Eq.([T 
and Ea.(|20p. we again find somewhat different scaling with the bias for the two failure prob 
abilities. For the Z failure probability the prefactor is 



V2~t 



■. exp 



2p 2 z mn 



e 4p z n _ i 



ozn -A ^~l 2 



2pzn\ 
TxZm ) 



1 y2ra 
8 n 



(e z -l) 1/2 ; (46) 



making the approximation e z — 1 w e z as before, and using m/n = \Jb = Ze z , this becomes 



2 Pz 

7T 



b -l/A Z -l/2 e Z'/S e Z/2 



2 Pz 

IT 



Z~ L e 



(47) 



For the X failure probability the prefactor is 



'27m 



exp 



2p 2 x mn 



e 4p x m _ i 



(e Apxm - 1) 



1/2 _ ( 2pz_ 

~TT~Z 



1/2 r 1 y2 n -i 



exp 



8 m 



e*-l 



(e* - 1) 



1/2 



(48) 



making the approximation e x — 1 ~ X as before, this becomes 
V 2 / v \ /o„ \ V 2 



"J Z CXP li^ 



We also note that, as in Eq. ([T"7j) . there are corrections to the factor e z / 8 of the form 

1 



exp -Z A [1 + 0(p z ) + 0{p z Z)] 



(50) 



In fact we can sum up these corrections to all orders in pzZ (still neglecting corrections 
suppressed by further powers olpz) to obtain the improved approximation 



e X p(^Z 2 (l-e-^)[l + 0(pz)}) 



(51) 



This correction can be significant for pZ = 0(1), but in that case our bound on the failure 
probability is rather loose anyway (see below), so including the correction is not so important. 
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Combining the estimated prefactors with Eq. (|45j) we obtain the formulas for the failure 
probabilities: 



BSZOptimalFailProb(pz, b) 
BSXOptimalFailProb{pz, b) 



1/2 



ZVb 
Z \ 8 PZ 



(^f\ 1/2 b-^(l-exp(-e- z ))^, (52) 



where Z = W(Vb). These asymptotic formulas apply if we fix the bias b at a large value and 
then allow pz to become sufficiently small. We see that the nonleading contribution to the 
log of the failure probability is approximately | Z 3 for Z-errors and approximately — \ In b for 
X errors. These nonleading terms are small compared to the leading term g^-Z 2 , which we 
minimized to obtain our estimates, provided pzZ R2 pzhxb <C 1- A numerical fit indicates 
that the leading correction to our asymptotic formula is approximately (.01)pz(lnb) s where 
i5«4. When this correction is small, we expect our estimate of the optimal failure probability 
to be reasonably tight. 

4 Comparison with numerics 

To check the accuracy of our formulas, we have also numerically determined the optimal 
values of m and n, and the corresponding optimal failure probability, for a variety of values 
of pz and b. 

Fig. [5] shows the numerically optimized failure probability and the asymptotic estimate 
Eq. (f23]) as a function of p = pz = px for the unbiased case (6 = 1). Our formula overestimates 
the exact result by less than 1% for p < .001 and by less than 10% for p < .01. 

BSOptimalFailProb( p) 



io- » 

10- !S 

io- 27 
io-' 6 



0.002 0.004 0.006 0.008 0.010 

Fig. 2. The optimal failure probability for the Bacon-Shor code (solid blue curve) and the estimate 
Kq. H23t (dashed red curve) for unbiased noise, as a function of p = pz = px- The curves nearly 
coincide. 

Fig. |3] shows the numerically optimized failure probability and the estimate Eq. (|52|) as a 
function of the bias b for pz — -01. Here we have actually plotted the total error probability 
assuming X and Z errors are independent; that is, 

BS Optimal F ail Prob = B S Optimal Z Fail Prob + BS Optimal XF ail Prob 

-BSOptimalZFailProb x B S Optimal X Fail Prob. (53) 
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The agreement is good for b > 10, though the discrepancy grows with increasing b, reaching 
about 10% for b = 500. 



BSOptimalFailProb( 0.01, b) 
II) 




Fig. 3. The optimal failure probability for the Bacon-Shor code (solid blue curve) and the estimate 
Eq J52H (dashed red curve) for pz = .01, as a function of the bias b = pz/px- The curves nearly 
coincide. 



Fig. 2] shows the numerically optimized failure probability and the estimate Eq. ([5^|) as 
a function of the bias b for pz = .03. Now our formula badly overestimates the failure 
probability for large 6, with the discrepancy reaching a factor of 13 for b = 5000. In this 
regime, the condition (.01)p z In 4 b -C 1 is not well satisfied ((.01)pz In 4 b w 1.6 for b = 5000). 



BSOptima!FailProb( 0.03, b 
10 




1000 2000 3000 4000 5000 

Fig. 4. The optimal failure probability for the Bacon-Shor code (solid blue curve) and the estimate 
Eq J52H (dashed red curve) for pz = .03, as a function of the bias b = pz/px- 



Fig.[S]shows the numerical and analytic results as a function of pz for b = 1000, and Fig. [5] 
shows the ratio of the two, illustrating that the agreement improves rapidly as p decreases. 
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BSOptimalFailProb( p. , 1000) 




0.02 0.04 0.06 0.08 0.10 

Fig. 5. The optimal failure probability for the Bacon-Shor code (solid blue curve) and the estimate 
Eq, l|52[l (dashed red curve) for b = 1000, as a function of the error probability pz- 



r( p- , 1000) 

4 - 



.1 - 



0.02 0.04 0.06 0.08 0.10 

Fig. 6. The ratio of the optimal failure probability to the estimate Eq, l|52j l for b = 1000, as a 
function of the error probability pz- The kinks in the plot occur because the optimal dimensions 
mxnof the Bacon-Shor code are integers which change discontinuously. 



5 2D Bacon-Shor recovery with syndrome measurement errors 

Measurement of the weight-two gauge operator XX can be executed by a simple circuit 
with four locations — e.g. preparation of an ancilla qubit in the X = 1 eigenstate |+), two 
successive controlled- A gates with ancilla qubit as control and data qubit as target, followed 
by measurement of the ancilla qubit in the X basis. A similar circuit measures ZZ, but with 
the controlled- A gates replaced by controlled- Z gates. 

A fault in the circuit that flips the ancilla qubit could cause the measurement outcome to 
be incorrectly recorded, and a faulty two-qubit gate could damage a qubit in the code block. 
However, a single fault in the circuit will not cause two errors in the code block that cannot 
be gauged away. A Z-type syndrome bit is obtained by computing the parity of n outcomes 
of XX gauge qubit measurements, and an error in any one of those n measurements could flip 
the bit, so that the probability of an error in a syndrome bit is roughly n times larger than 
the probability of error in each gauge qubit measurement. Similarly, an A-type syndrome bit 
is obtained by computing the parity of m outcomes of ZZ gauge qubit measurements, so the 
probability of an error in a syndrome bit is roughly m times the probability of error in each 
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gauge qubit measurement. 

Though measurement and data errors may actually be correlated, these considerations 
motivate an idealized model for Bacon-Shor error recovery with noisy measurements. For the 
Z-type syndrome we consider a repetition code of length m, where both the probability of 
error per data bit and the probability of error per syndrome bit is pz = OddProb(pz,n), 
and pz is the fixed physical error rate for a Z-type error. Likewise, for the X-type syndrome 
we consider a repetition code of length n, where both the probability of error per data bit 
and the probability of error per syndrome bit is px = OddProb(px , tn), and px is the fixed 
physical error rate for an X-type error. 

We will study this model for both unbiased noise and highly biased noise with b = 
Pz/px ^ 1- Justifying highly biased noise models raises subtle issues — why should the 
nontrivial quantum gates used in the syndrome measurement circuit strongly favor Z errors 
over X errors [TO] [9]? We will not address these issues here; rather we shall be satisfied to 
say that a model in which both the data qubit errors and the syndrome bit errors are highly 
biased is mathematically natural and worthy of investigation. Ref. 9; contains a much more 
complete discussion of fault-tolerant error correction for Bacon-Shor codes. 

A nearly optimal recovery scheme for a repetition code with syndrome errors was described 
in [TTJ [T^] . To obtain more reliable syndrome information we follow the history of the syn- 
drome through many measurement cycles, and assume that the probabilities for both data 
errors and syndrome errors are given by pz, Px in each round of syndrome measurement. 
The syndrome history can be represented on a two-dimensional square lattice in spacetime, 
where each horizontal row of squares records the results from one cycle of (possibly noisy) 
syndrome measurement and each square represents a data bit — a marked vertically oriented 
link in the row indicates that the two data bits sharing that link were found to have opposite 
parity. We identify all the boundary points of this chain of nontrivial syndrome bits, and 
apply the Edmonds matching algorithm to find the minimum-weight chain in spacetime with 
these boundary points (relative to the boundary of the sample). This minimum- weight chain 
identifies the most likely error history compatible with the observed syndrome history, where 
vertical edges in the chain are hypothetical syndrome measurement errors and horizontal 
edges are hypothetical data errors. The actual error chain combined with the hypothetical 
chain comprises a closed chain relative to the boundary, and a segment of this chain stretch- 
ing across the sample, connecting it's left and right boundary, signifies a logical error. If we 
consider the repetition code on a circle rather than an open one-dimensional lattice, so there 
is no sample boundary, then a closed chain wrapping around the cylinder indicates a logical 
error. 

This scheme has been studied previously [Til HI] in the case where the probability p of a 
data error or syndrome error is a constant independent of the length m of the repetition code; 
in that case recovery is successful with a probability approaching 1 as m — > oo for p < 10.3% 
using the matching algorithm, and for p < 11.0% using the optimal algorithm. Now we wish 
to reconsider the efficacy of the matching algorithm in the case where the error rate scales 
nontrivially (roughly linearly) with m. 

Based on our earlier results for the case of ideal syndrome measurement, we anticipate 
that if m is chosen optimally then the probability per unit time of an encoded error in the 
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Bacon-Shor code has the form 

ln[BS Optimal F ail Rate(p z , b)] « -A(b)/p z + B{b) lnp z + C(b) + ■ ■ ■ , (54) 

but with different functions of the noise bias A(b), B(b), C(b) than in the case of ideal mea- 
surements. We could try to estimate these functions using a Monte Carlo method, in which we 
generate sample error histories, infer the syndrome, find the minimum weight matching of the 
syndrome's boundary points, and then determine whether a logical error has occurred. This 
method is difficult to carry out, however, because for small pz the optimal failure probability 
is quite small, and we need to generate many samples to estimate it with reasonable statistical 
accuracy. Instead, we will use an analytic argument to obtain a rather loose upper bound on 
A (which dominates the scaling of the failure probability when pz is small) for both the case 
of unbiased noise (b = p z /pz = 1), and the case of highly biased noise (6 = pz/px 3> !)• 

Since we expect the value of A to be the same for a planar Bacon-Shor code as for a code 
defined on a torus, we will consider the case of a torus for convenience. Thus we consider the 
syndrome history (for either the X or Z errors) on a cylinder, closed in the spatial direction 
but open in the temporal direction. If there is a logical error, then the combination of the 
actual error chain and the hypothetical error chain must contain a self-avoiding cycle that 
wraps once around the cylinder, where at least half of the edges in this cycle have actual 
errors (otherwise there would a lower- weight choice for the hypothetical chain). If we fix a 
cycle of length r, then the probability that a particular set of s edges in the cycle have actual 
errors, while the remaining r — s edge do not, is 

f(i-pys = [p(i-p)Y/2[p/(i-p)y-r/2 < [p(i-p)f/ 2 (55) 

where the inequality holds provided p < 1 — p and s > r/2. The number PJ of ways to place 
m errors on the cycle is bounded above by 2 r , therefore, the probability of failure arising from 
this particular length r cycle is bounded above by [4p(l — p)} r / 2 . 

Considering the Z error syndrome history for definiteness, we now let SAC(r, m) denote 
the number of distinct self-avoiding cycles of length r that wrap around the cylinder of cir- 
cumference m. Here we are not counting the freedom to translate the cycle in either space 
or time. Taking this freedom into account, we obtain an upper bound on the probability of 
error per unit time 

oo 

BSZFailRate{p z ,m,n) <m^] SAC(r, m)[4p z (l ~ pz)Y /2 - (56) 

r— m 

The factor of m in front arises from the m possible spatial translations of the cycle, and 
a factor of time T arising from time translations has been divided out to obtain a failure 
rate per unit time. Note that the cycle length is in principle unbounded, as it could extend 
indefinitely in the time direction. 

The number SAC (r, m) of self-avoiding cycles of length r on the cylinder is bounded above 
by the number SAW{r) of self-avoiding open walks with a specified starting point, for which 
an upper bound is known of the form [131 114] 



SAW(r) < 7 rV, 



(57) 
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where fj, w 2.6381585, f3 = 11/32, and 7 w 1.17704. Plugging the upper bound on SAW(r) 
into our expression for BSZFailRate, we find 

00 

BSZFailRate(p z ,m,n) < 7 m^ r f3 fi r [4p z (l - Pz)T /2 

r—m 

s—m ^ ' 

00 

< im f i+1 {4n 2 p z (l - p z )] m ' 2 Y, e tis ' m [^ 2 p z {l - Pz)] s/2 



s=0 

1 - V 4 M e 2 p/ m p z {l - pz) 

We will obtain our upper bound on the failure rate by choosing m and n such that 4/i 2 e 2 ^^ m p z (l— 
p z ) is strictly less than one; therefore we can bound the last factor by a constant, obtaining 

\n[BSZFailRate(p z , m, n)] < y ln[4^ 2 p z (l -pz)}+ 0(ln m), (59) 
and by similar reasoning 

\n[BSXFailRate(p x , m, n)} < | ln[4 M 2 p x (l - p x )] + O(lnn). (60) 
5.1 Unbiased noise 

In the case of unbiased noise, where pz = px = P, wc choose n = m and recall that p = 
\ (1 - (1 - 2p) m ) sa \ (1 - e- 2 P m ), so that 

4p(l-p) w 1 -e- 4pm . (61) 

Thus we conclude that the failure rate for either Z-type or X-typc errors is 

ln[BSFailRate(p,m)} < j In [^ 2 (l - e - 4pm )] + O(lnra). (62) 

To find the value of m that minimizes this upper bound (for asymptotically large m), we 
solve 



= \ (m [M 2 (1 - e- 4 -)] + + 0(l/m), 



(63) 



for fi = 2.6381585, finding pm w 0.0139682, so that 



00679079 

\n[B S Optimal F ail Rate(p)] < + O(lnp). (64) 

P 

Thus we derive a lower bound 

A > .00679079, (65) 

which is about an order of magnitude smaller than the value of A found for the case of ideal 
syndrome measurement. 
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Even more simply, we may observe that p < pm and p(l — p) < p to obtain 

]n[BSFailRate(p,m)} < y ln(4/x 2 pm) + O(lnm), (66) 

which is optimized by choosing 4p 2 pm = 1/e, so that 

1 , .00660714 

8/i 2 ep 



\n[BSOptimalFailRate{p)] < —^-^ + O(lnp) = + O(lnp), (67) 



which provides nearly as good a lower bound on A. 
5.2 Highly biased noise 

In the case of biased noise, again using the approximations pz(l~pz) < Pzn andpx(l— Px) < 
pxTn, the leading contributions to the Z and X failure rates are 

Tn • I 1 
\n\BSZFailRate(pz,ni,n)] w — Infii^pzn) — — 5 XlnZ = — 5 — bX\nZ, 

ln[BSXFaiZi?afc(px, m, n)] ~ ^ln(4/i 2 p x m) = — - — ZlnX, (68) 

Z ofl pz 

where 

X = 4/i 2 p x m, Z = 4/i 2 p z n, b=—. (69) 

(Actually, the approximation pz(l — pz) * Pz^ is not necessarily so accurate when m and n 
are chosen optimally, but it suffices for deriving an upper bound on the failure rate; if we do 
not make this approximation the analysis becomes more complicated and the upper bound 
improves by only about 8%.) Equating the upper bounds on the X and Z failure rates, we 
are to minimize 

F{X,Z) = Z\nX (70) 

subject to the constraint 

ZlnX = bXkiZ. (71) 

Introducing a Lagrange multiplier A to do the constrained minimization, we obtain the equa- 
tions 

(l + \)kiX = \b—, 
(1 + A)4 = A&lnZ > ( 72 ) 



which imply 



lnX = i i_. (73) 



Using the constraint Eq. (jTTj) , we see that the aspect ratio of the optimal code is 
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When the bias b is very large, the optimal value of Z (that is, when m and n are chosen 
to optimize the approximate expression Eq. (|68|) . an upper bound on the actual failure rate) 
approaches one from below, while 4/i 2 e 2 ^/ m pz(l — pz) remains strictly less than one for m 
large, so we can still bound the last factor in Eq. (f58f by a constant. To estimate the failure 
probability, we first determine Z by solving Eq.([7Tj) in the form 

Z = b In 2 Z exp(l/lnZ). (75) 

Then we find X by substituting this value of Z back into Eq. (|7T|) . finding 

b \n 2 X 1 , / 1 \ l~b 



=\T^- (76) 



z x \Vx) V 4Z 

Recalling that the Lambert W function W(z) is defined by We w = z, we conclude that 



- ' x - 

Thus we find that the value of F at its minimum is 

F mia = Z]nX^-2Zw(J^J , (78) 
and we obtain upper bounds on the leading contributions to the failure rates: 



\n[BSZFailRate{pz,m,n)} w ln[BSXFailRate(p x , m, n)] - W J , (79) 

where Z(b) is a slowly varying monotonic function, ranging from 1/e to 1 as b increases from 
1 to infinity. Comparing to the upper bound for the case of ideal syndrome measurement we 
find 

a slowly increasing function of b. The optimal aspect ratio is 



in contrast to the case where the syndrome is perfect, this aspect ratio grows only polyloga- 
rithmically with the bias b. 

For example, by solving Eq.^), we find Z = .3679, .7046, .8428, .8994 for 6 = 1, 10 2 , 
10 4 , 10 6 and corresponding values of -F min = .3679(= 1/e), 2.01255, 4.92966, 8.48426. For 
b = 10 2 , 10 4 , 10 6 , our values of of A/ A are 5.231, 11.50, and 18.29 respectively. 
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6 Higher-dimensional Bacon-Shor codes 

The two-dimensional Bacon-Shor code can be extended to higher dimensions. Bacon 1 
discussed one such extension to three dimensions. In Bacon's code, defined on an L x L x L 
cubic lattice (for L odd) , the logical operators X and Z are both weight L 2 operators defined 
on orthogonal planes, and a length L repetition code protects against both X and Z errors. 
Here we will briefly discuss a different way to extend the code to three (or more) dimensions, 
in which the logical operator Z resides on a horizontal plane and the logical operator X 
resides on a vertical line. This alternative formulation is less symmetric but in some ways 
more natural; it may be more effective than Bacon's code if the noise is highly biased, and it 
is more robust against errors in the measurement of the Z-type error syndrome 

We consider an m x n x k cubic lattice, and choose the gauge algebra to be generated by 
two-qubit operators XX acting on pairs of neighboring qubits in the same horizontal plane, 
and by two-qubit operators ZZ acting on pairs of qubits in the same vertical column. There 
is one protected qubit; the weight-mn operator Z acts on all qubits in a horizontal plane and 
the weight-fc operator X acts on all the qubits in a vertical column. The stabilizer group has 
ran — 1 independent X-type generators, each a product of 2k Xs acting on all the qubits in a 
pair of neighboring vertical columns, and k — 1 independent Z-type generators, each a product 
of 2mn Zs acting on all the qubits in a pair of neighboring horizontal planes. As in the 2D 
code, the value of a check operator can be found by measuring weight-two gauge operators 
and combining the results. For example, the weight-2mn product of Zs in two neighboring 
planes is obtained by measuring ran gauge operators ZZ, each acting on two qubits in the 
same column, and multiplying the results together. 

As in the 2D case, the error recovery procedure is clarified if we adopt a suitable gauge. 
An even number of Z errors in a column can be gauged away, and if there are an odd number 
of Z errors in a column, there is a gauge such that all errors but one are removed, and the 
single error lies in the uppermost horizontal plane. Similarly, an even number of X errors in 
a plane can be gauged away, and if there are an odd number of X errors in a plane, all but 
one can be removed, with the single error lying in the northwest corner of the plane. The 
check operators, whose value can be inferred from measurements of many weight- two gauge 
operators, check whether qubits in the northwest corner of neighboring planes agree in the Z 
basis, and whether qubits in the top plane of neighboring columns agree in the X basis. 

The planar repetition code corrects (run — l)/2 Z errors in the plane, and the columnar 
repetition code corrects (k — l)/2 X errors in the column. Aside from its greater length (as- 
suming mn > fc), the planar code has another advantage over the columnar code — measuring 
all local X-type gauge operators determines the error syndrome redundantly, which improves 
the robustness against syndrome measurement errors. From the gauge qubit measurements, 
we determine the value of X® 2k for any two neighboring columns; there are m(n — 1) such 
pairs of "east-west" neighbors as well as n(m — 1) pairs of "north-south" neighbors, for a total 
of 2mn — to — n check bits, or which only mn — 1 are independent. 

Recovery from Z errors succeeds if the number of vertical columns with an odd number of 
Z errors is no more than (mn — l)/2. and recovery from X errors succeeds if the number of 
horizontal planes with an odd number of X errors is no more than (fc — l)/2. Setting m = n, 
therefore, the probability of a logical Z error in an independent error model (assuming perfect 
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syndrome measurements) is 

BS3DZFailProb(p z , m, k) = RepFailProb(OddProb(p z , k), m 2 ) (82) 
where pz is the Z error rate, and the probability of a logical X error is 

BS3DXFailProb(p x , m, k) = RepFailProb(OddProb(p x ,m 2 ), k) (83) 

where px is the X error rate. These are the same functions we encountered in our analysis 
of the 2D code, except with m replaced by m 2 and n by k. 

In four dimensions, denote the four directions by x,y,z,w. To define the Bacon-Shor 
code, choose the gauge algebra to be generated by XX acting on neighboring qubits in each 
xy-plane (with z and w fixed) and by ZZ acting on neighboring qubits in each zw-plane 
(with x and y fixed). Note that each xy-plane intersects exactly once with each zw-plane. 
The logical Z is the product of all Zs in an zy-plane, and the logical X is the product of all 
Xs in a zw-plane. For a suitable gauge choice, all X errors can be moved to the ^w-plane 
with x = y = and all Z errors can be moved to the xy-plane with z = w = 0. Hence for this 
four-dimensional Bacon-Shor code a two-dimensional repetition code protects against logical 
Z errors and another two-dimensional repetition code protects against logical X errors. If we 
choose the dimensions of the hypercubic lattice to be m x m x n x n, then the probability of 
a logical Z error (assuming perfect syndrome measurement) in an independent noise model is 

BS4DZFailProb{p z , m, n) = RepFailProb(OddProb(p z , n 2 ), m 2 ) (84) 

where pz is the physical Z error probability per qubit, and the probability of a logical X 
error is 

BS4DXFailProb(p x ,m, n) = RepFailProb(OddProb(p x ,m 2 ), n 2 ), (85) 

where px is the physical X error probability per qubit. These are the same functions as for the 
2D code, but with m replaced by m 2 and n by n 2 . The advantage of the four-dimensional code 
over the two-dimensional or three-dimensional code is that measuring the weight-two local 
gauge operators provides redundant syndrome information (and hence improved robustness 
against syndrome measurement errors) for both the Z and X errors. 

7 Conclusions 

We have studied the performance of Bacon-Shor codes against both unbiased and biased 
noise, concluding that the codes provide excellent protection against logical errors if the error 
probability per qubit is less than a few tenths of a percent in the unbiased case, and if the 
dephasing error probability is less than a few percent in the case of highly biased noise, 
assuming the syndrome is measured perfectly. 

Using quantum codes to protect quantum computers from noise raises many thorny prob- 
lems 015]; in particular, we need to build a set of gadgets that reliably execute logical gates 
acting on the code space, and in the case of biased noise we need to employ a physical gate 
set compatible with the bias [TU1IH]- We have not discussed these issues here. Instead, our 
goal has been to gain a better understanding of the properties of the Bacon-Shor code family, 
without getting bogged down in detailed constructions of fault-tolerant protocols. We expect, 
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though, that our analytic formulas for the logical failure probabilities of optimal codes will 
provide helpful guidance for the construction of fault-tolerant schemes. 
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